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SEARCHING FOR DIOPHANTINE QUINTUPLES 


MIHAI CIPU AND TIM TRUDGIAN 


Abstract. We consider Diophantine quintuples {a, b, c, d, e}. These are sets of distinct 
positive integers, the product of any two elements of which is one less than a perfect 
square. It is conjectured that there are no Diophantine quintuples; we improve on current 
estimates to show that there are at most 1.18 • 10^^ Diophantine quintuples. 


1. Introduction 

Define a Diophantine m-tnple as a set of m integers {oi,..., am} with oi < 02 < ... < 
Qm, snch that OjOj + 1 is a perfect sqnare for all 1 < z < j < m. Thronghont the rest of 
this article we frequently refer to m-tuples, and not to Diophantine m-tuples. 

It is conjectured that there are no quintuples — see [21 [13]. Successive authors (see, 
e.g., Table 1 in [16]) have reduced the bound on the possible number of quintuples. The 
best such published bound is 2.3 • 10^® by Trudgian [16]. The purpose of this paper is to 
improve on this in the following theorem. 

Theorem 1. There are at most 1.18 ■ 10^’^ Diophantine quintuples. 

In ^we collect some ancillary results that aid the computational search for quintuples. 
In l|3]we obtain bounds on the relative sizes of elements in a quintuple. We use this in ^ 
with results on linear forms of logarithms to obtain upper bounds on the second-largest 
element in a quintuple. In ^we examine some number-theoretic sums, which enable us to 
bound the total number of quintuples. We present two new arguments in ^that enable 
us to make a further saving, and ultimately to prove Theorem [T] 

We are grateful to Adrian Dudek who suggested the asymptotic form in (jHI). 

2. Discards 

It is known that every triple {a, b, c} can be extended to a quadruple of a certain form. 
This is dubbed the ‘regular’ quadruple and is denoted as {a,b,c,d+}. If a double or a 
triple cannot be extended to a non-regular quadruple, then it cannot be extended to a 
quintuple. We call such doubles or triples discards. The doubles {k, k + 2} [10] (see also 
0) are discards for k > 1. For an extensive list of discards, one may see [161 §2-1]. The 
following result allows us to recognise many discards. 
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Lemma 2.1. Let {a, b, c, d} be a Diophantine quadruple with a<b<c<d^<d. 

• Ifb<2a then b > 21000. 

• If 2a<b< 12 a then b > 130000. 

• If b > 12 a then b > 4001. 

Proof. The only difference between this lemma and Lemma 3.4 in [5] is the exclusion of 
the value b = 4001 in the last case. Indeed, a pair {a, 4001} with 12a < 4001 cannot 
be extended because the equation 4001a + 1 = has unique integer solution r < 4001, 
namely r = 4000, which entails a = 3999. □ 

Lemma 2.2. ([SI Theorem 1.1]) Let {a, 6, c, d, e} he a quintuple with a<b<c<d<e 
and put g = gcd(a, 6). Then b > Sag. If moreover c > a + b + 2\/ab + 1 then b > 
max{24 ag, 2 a^^‘^g‘^}. 

Lemma 2.3. ([5l Theorem 1.3]) Let {a,b,c,d,e} be a quintuple with a<b<c<d<e 
and c = a + b + 2\/ab + 1. Then b < a? and gcd(6, c) = 1. In particular, at least one of 
a, b is odd. 

Examination of the relative size of entries in a quintuple has the following outcome. 

Lemma 2.4. Any quintuple {a,b,c,d,e} with a<b<c<d<e must be of one of the 
types listed below: 

(A) 4a < 6 and Aab + b + a < c < 6^/^, 

(B) 4a < 6 and c = a + b + 2\/ab + 1, 

(C) 4a < 6 and c > 6^/^, 

(D) 6 < 4a and c = a + b + 2\/ab + 1. 

Proof. The claim follows from the previously cited results having in view 0 Lemma 4.2] 
and Subsection 2.2 in |1]. □ 

3. Exploiting the gonnegtion with Pellian equations 

The entries in a quadruple are severely restricted in that they appear as coefficients of 
three generalized Pell equations that must have at least one common solution in positive 
integers. Each component of such a solution is obtained as a common term of two second- 
order linearly recurrent sequences, giving rise to relations of the type z = Vm = Wn for 
some positive integers m and n. A key ingredient in the study of Diophantine sets is a 
relationship between the parameters m, n, and the values in the set in question. 

Our next result is of this kind. It improves on several versions already in the literature 
— see, e.g., jHHSlIIZ]. 

Proposition 3.1. Let {A, B,C, D} be a quadruple with A < B < C < D for which 
V 2 m = W 2 n has u solutiou with 2n > m > n > 2, m > S. Suppose that A > Aq, B > Bq, 
C > Co, B > pA for some positive integers Aq, Bq, Cq, and a real number p > 1. Then 

m > aB-^Cc^P^ 
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where a is any real number satisfying both inequalities 

(1) ^Cq < 4, 

(2) + (4Eo(A + + 2(A + a < ABq, 

with A = (Aq + + 1)“^/^. 

Moreover, if > f3B for some positive real numbers jS and r then 

m > 

Proof. We assume that m < aB~^B(j^P aim at establishing a contradiction if a is 
too small. We start from the congruence (see, e.g., [H Lemma 4]) 

(3) eAimf + Sm = eBn^ + Tn (mod 4(7) for some £ = ±1, 
where S = \JAC + 1 and T = \JBC + 1. Since 

{Am"^ — Bn^\ < maxIAm^, Bn^} < Bmf < o? C 

and 

|S'm — Tn\< max{S'm,Tn} < Tm < aB~^^‘^C^^‘^VBC + 1 

< ai?-V2c'V2(5i/2c'i/2 + 15-1/2^-72) <a{l + lBf^Cf^)C, 

then, if a satisfies ([1]), the congruence ([3]) becomes the equality Am^ — Bn^ = e{Tn — Sm). 
Multiplication by Tn + Sm followed by rearrangements results in the equality 

(4) (Tn^ — Am'^) ((7 + e{Tn + S'm)) = m^ — n^. 

Note that Bn^ = Am^ entails m^ = n^, so that A = B\ a. contradiction. Hence, for m = n 
one necessarily has (7 = Tn + S'm, while for m > n one finds that Tn^ — Hm^ divides 
the positive integer m^ — n^, so that m^ — n^ > |Hm^ — Tn^|. This gives the following 
inequality 

m^ T + 1 

— > - . 

n2 A +1 

Having in view the lower bounds for A and B, we obtain 

m^ pA + 1 pAq + 1 1 

n2 A + 1 ~ Aq + 1 A2 

From dlj), m < 2n, and the definitions of S and T, we conclude that 

(7 < Tn + Sm + m^ — n^ < \m\/ BC + 1 + mV AC + 1 + | m^ 

< I a^B-^C + aB-^/^C^'^ [xVBC + 1 + Vp~^BC + l) 

< I + aC- (A(l + \B-^C-^) + p-^/\l + ipT-^C'-')) 

< I a^Bf^C + aC (A(l + \Bf^Cf^) + + \pBf^Cf^)) . 
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The last expression is at most C if a satisfies the inequality ([2]), whence the hrst inequality 
in the conclusion of our proposition. The second one is readily obtained from what we 
have just proved and the hypothesis > f3B. □ 

Lemma 3.1. If {a,b,c,d,e} is a quintuple with a<b<c<d<e then the following 
bounds for m hold: 

(A) m > 3.3022di/^ (B) m > 1.5002^2/7^ (C) m > 2.0604^3/10^ 

(D) m > 1.0080^3/3. 

Proof. This is an application of the result just proved for {A, B, C) = (a, b, d) in cases 
(A)-(C) and for (A, 5, (T) = {a,c,d) in the remaining case. We use Proposition 13.11 with 
carefully chosen values for parameters in ranges suggested by Lemmas I2.1H2.3I In [12] 
it was shown that the hrst hypothesis required in Proposition 13.11 is satished. It is also 
known that one has d > Aabc + a + b + c (see, for instance, the proof of Lemma 6 in [S]). 

In case (A) Lemma 12.21 hints to consider separately values of a less than 144 since then 
one has B = b > max{24a, 20 ^/ 2 } = 24a = 24A. A short computer search hnds that 
Bq = 4095 and ^ > 146 =: p. Clearly we must put Aq = 1. From 

C = d > Aabc + a + b + c> {Aab + l)(4a6 + a + b) > (16a^ + 4a)6^ 

it follows that r = 1/2, (3 = (16Aq + 4Ao)3/2^ Cq > 3.35 • 10® are admissible choices. Both 
inequalities ([1]) and ([2]) are satished by a = 1.56155. 

Still in case (A), when a > 144 one puts Aq = 144, Bq = 4002 (by Lemma [2All . p = 24 
(see Lemma 12/21) . r = 1/2, fl = (16Ag+ 4Ao)3/2^ whence Cq > 5.32-1032 and a = 1.56155. 

Having in view Lemma 12.11 in case (B) we hrst examine the subcase 4a < 5 < 12a. 
Then Bq = 130001, which implies Aq = 10834 and p = 4. From 

c> 6(1 + 12-3 ^ 2 .12-3/2) = (1 + 12-3/2)2.0 
and a® > 6 it follows that 

C = d> {Aab + l)(a + 6 + 2r) > 4(1 + 12 - 3 / 2 ) 2 a 62 > 4(1 + 

so that r = 3/7, /3 = (2 + 3 - 3 / 2 ) 6 / 7 ^ Cq = 5.68 ■ 1032. For these choices it is readily 
obtained that a = 0.9999 is permissible. 

The other possibility in case (B) is to have b > 12a. Convenient values of parameters 
are p = 12, Aq = 16 (from > b > 4000), Bq = 4002, r = 3/7, (3 = 2^l\ Cq = 1.01 ■ lOT 
for which the same value a = 0.9999 works. 

Case (C) is similar to case (A). Now, for a < 143 we see that we can take Aq = 1, 
Bq = 4004, p = 28. As 

C > 4a6c > 4a6^/2 > 4.05 ■ 10^ =: Co, 

we further get r = 2/5, f3 = 42 /^^ whence again a = 1.5615. In the complementary 
subcase a > 144, admissible values are Aq = 144, Bq = 4002, p = 24, r = 2/5, /3 = 5762/^, 
Co = 5.83 • 1033. Plugging these specializations into Proposition 13.11 we obtain the same 
value for a. 
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Finally, in case (D) we have A = a < 6/3, B = c = a + b + 2\/ab + 1 > (1 + 

B < a + b + 2^?)~^b{b — 1) + 1 < (1 + 3“^/^)^6, and 

C = d> 4abc > b^c > (1 + 

Therefore, p = (1 + 3^/^)^, r = 1/3, and /5 = (1 + 3“^/^)“"^/^. From 130001 < 6 < 4a, we 
have Aq = 32501, whence Bq > 292504, and Cq > 4.04 ■ 10^^. From ([1]) and ([2]) we obtain 
a = 1.3660. □ 

For future reference, the values used in the previous proof are given in Table [TJ 


Type 

Aq 

Bo 


P 

f3 

T 

(AI) 

1 

4095 

3.35 • 10« 

146 

201/2 

1/2 

(All) 

144 

4002 

5.32 • 10^2 

24 

24 • 5771/2 

1/2 

(BI) 

10834 

130001 

5.68 • 10^2 

4 

(2 + 3-1/2)6/7 

3/7 

(BII) 

16 

4002 

1.01 ■ 10^ 

12 

26/7 

3/7 

(Cl) 

1 

4004 

4.05 ■ 10^ 

28 

42/5 

2/5 

(CII) 

144 

4002 

5.83-1011 

24 

5762/5 

2/5 

(D) 

32501 

292504 

4.04 • lOii' 

(1 + 31/2)2 

(1 + 3-1/2)-4/3 

1/3 


Table 1. Parameter values for various types of Diophantine quintuples. 


The values of a, and hence the bounds on m in Lemma [3.11 rely on the computational 
bounds in Lemma 12.11 While it is tempting to extend these computations, such an 
extension would have almost no effect on the values of a. Consider, for example, case 
(A): sending Bq, Cq to inhnity in ([T]) gives + a < 4. Therefore the optimal value of a 
is 1.5615528..., whereas we have a = 1.56155. Likewise, in case (D) the optimal value is 
i(l + \/3) = 1.366025..., whereas we have 1.3660. It seems that a new idea is needed to 
improve substantially on the lower bounds on m. 


4. Employing linear forms in the logarithm 


The lower bounds for the index m given in the previous section can be complemented 
by inequalities derived from upper bounds for linear forms of logarithms of algebraic 
numbers. To this end, we apply the best result of which we are aware. 


Theorem 4.1 (Aleksentsev). Let A he a linear form in logarithms of n multiplicatively 
independent totally real algebraic numbers Oi,... with rational coefficients bi,... ,bn- 
Let h{aj) denote the absolute logarithmic height of aj for 1 < j < n. Let d be the degree 
of the number field /C = Q(ai,..., an), and let Aj = ma.x{dh{aj), \ loga^l, 1). Finally, let 


(5) 

Then 


E 


max ( max 

\l<i,j<n 


N , N\ 

41, ^ A, / 



log|A| > —5.3n + l)"'’''^(n + 8 )^(n + 5)(31.44)"'(i^(logE)Ai • • • A„log(3n(i). 
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We have used the hrst displayed equation on [U p. 2] to dehne i? in (|5]): this makes our 
application easier. We apply Theorem 14.11 for d = 4, n = 3 and to 

A = j log ai — k log 0^2 + log tta, 


with 

\/A ± 

where the signs coincide. 

For our purposes we do not need the exact values of Aj and E as dehned in Theorem 
su decent estimates will suffice. To hnd these estimates we proceed as follows, keeping 
the notation and hypotheses of Proposition 13.11 and supposing additionally that C < Ci 
for a certain integer Ci. 

We begin by noting that one has 

21ogai < log(4AC' + 4) < log(4p“^(i? — 1)C + 4) < log(4p“^i?C') 

< log(4p-V-'C'+") 

provided that pA <5 — 1. This clearly follows from pA < B when p is integer, as in 
cases (A)-(C). In case (D) we have 6 > 3a + 1, so that (cf. the proof of Lemma ITT]) 

B = c = a + h + 2\/ ah + 1 > 1 + (1 + = 1 + pA. 

In each of the cases (A)-(D) we have f3p > 4, whence 

Ai < 5 (i(/?,p,r,Ci) logC, 


«! = S ' + Vac, a 2 = t + Vbc, 03 



with 


gi{(3,p, T, Cl) 


log4-log(/?p) 

1 ^ -i-7^- 

log Cl 


We readily obtain the following lower bound on Ai 


with 


Ai > log(4AC) > g 2 {Ao,Ci)\ogC, 

fA ^^log4 + logAo 

P2(Ao,Ci) := 1 + 


log Cl 

Similar relations hold for A 2 , namely 

2 log 02 < log(45C + 4) < \og{4:l3~^C^^'^ + 4), 
which implies the upper bound 


A 2 < g3(/3,T,e) logC, 


93(/^,C e) 


log4 + log(/3 ^ + e ^ 

1 + r H-^- 

log e 


where 













SEARCHING FOR DIOPHANTINE QUINTUPLES 


7 


and e = Co in the cases (B), (Cl), and (D) (when /? < 4) and e = Ci in the remaining 
cases (A) and (CII). An easily-derived lower bonnd for A 2 is 

^2 > fi' 4 (-Bo, Cl) logC, 


with 


The ineqnalities 


94(^0, Cl) 


log 4 +log Eq 
log Cl 


y/B \/C + \/A^\/B \/C + \/A^^ VB \/C-\/A^^ 

-^fB ^ 

are obvions. The modulus of the fourth algebraic conjugate of 03 is also greater than 1 
precisely when \/C(\/S — \/A) > 2\/ AB. This inequality holds whenever 

( 6 ) pBl-^{p^/^ - 1 ) 2 " > 22 ". 


It is easy to check that (jb]) is satished in each of the cases (A)-(D). One now obtains 

^ 3 . 4 h(a. 3 ) ^ log . 

where g is the content of the polynomial A2(C —i?)2X^-1-4^2i?(C — i?)X^-|-2Ai?(3Ai? — 
AC — BC — C2)X2 -|- 4Ai?2(C — A)X -|- B‘^{C — A)2. Since g is at most the smallest of 
the coefficients, which is 4A2i?(C — B), one has 

Note that B(C — A) < /5“^C^+" readily implies 


-43 < g5(/3,T,f)logC, 


with 


95{P,T,f) := 2 + 2 r- 


2log/3 

log/ 


and / = Cl if /5 > 1 and / = Cq if /5 < 1. A lower bound for A 3 is obtained with the help 
of the inequalities 


A 3 > log 


B{C -Af 
4A2(C-5) 


> log 


/jp2c'^-"(i - AoCr ^)2 
4(l-pAoCri) 


which entail 
where 


-43 > 5 ' 6 (/^, P, r, Ao, Cl) log C, 


ge{P,p,T,Ao,Ci) := 1 - t + 


log(|/jp 2 ) ^ 2 log(l - ApCr^) - log(l - 4Cr^) 
log Cl 


On noting that for all relevant values of parameters one has 5 ' 2 (Ao,Ci) < g 4 {Bo,Ci) 
and using the inequality g 2 > 9e (which follows, for Ci > 10 ^ 2 ^ from 16C2"(1 —/?“^) > j3p‘^ 
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if /3 > 1 and from 16(1 + 3 in case (D)) as well as the known 

relation k < j, we find that we may take 


E < 


2j 


96{/3,P, r, y4o,Ci)logC'o' 
Hence, Theorem 14 .1 1 yields the following corollary. 

Corollary 4.2. 

-logA < 1.5013 ■ 10 ^^ 5 ( 35 ( 5(2 logQ;i)(log^C') log 


2j 


96 log Co 


Corollary 14.21 bounds A from below; we can bound A from above using Eq. (4.1) in [12], 
which states that 

0 < A < -ACaf^. 

3 ^ 

Comparison with Corollary 14.21 gives the main result of this section. 

Proposition 4.3. 


j < 1.50131 ■ 10 535-5 (log C) log 




^56 log Co/ ■ 

Set j = 2m in Proposition 14.31 and use Lemma [3.II with the values given in Tabled] and 
Cl = 10 ^^-^®® in all cases, as per jH Theorem 1.2]. We thus get a new upper bound on d 
that we take as Ci in a new iteration of this procedure. Slightly better bounds result by 
taking much higher Co (just below the value for Ci considered in the same iteration). We 
record our computations in the following theorem. 

Theorem 2. If {a, b, c, d, e} is a quintuple with a<b<c<d<e then the following 
bounds for d hold: 

(A) d < 1067-859 < 7.228 ■ 1067, (B) d < 1060-057 < 1.141 . 1060 ^ 

(C) d < 1066-528 < 3.373 . iq56^ P) ^ ^ 1Q51.416 ^ 2.603 ■ lO^T 

We close this section with a remark concerning the size of the smallest entry in a quin¬ 
tuple. Although it has no immediate bearing on the next section, further improvements 
on d should enable future researchers to enumerate all possible triples. Recording the 
maximal size of a should aid this goal. 

Proposition 4.4. The only quintuples that could arise from case (A) are those in which 
a < 7.4 • 107. 

Proof. The triples in case (A) must satisfy 56/2 > c > 4a6 + 6 + a, so that, in particular 
a < Some quick computations give that for Aq = 7.4-10^ one obtains d < O.l-lO^^. 

From d > Aabc > IQafb'^ > (lOa^)^ it then follows a < 7.29 • 107, a contradiction. □ 











SEARCHING FOR DIOPHANTINE QUINTUPLES 


9 


5. Bounding the total number of quintuples 


In this section we combine the methods of |1] and [16] in bounding certain arithmetical 
sums. Let 


E{x) = 

n=l 


F(x) = 

n=l 


2^iu(n) 

n 


^ ocj(2n—1' 

G(^) = y H-—: 

^ ^ ^ 2 n - 1 

n=l 


where all the sums are defined for real values of x > 1. Bounds for these sums can be 
used as in the proof of Lemma 4.4 in |1] to prove that 


N 


(7) -1) <‘2NG 


n=2 


H + l 


+ NG 


77 + 4' 


+ NG 


77 + 2' 


+ NG 


77 


where duin) counts the number of divisors of n that do not exceed 77. The function 
dni^^ — 1 ) arises naturally when considering the number of doubles {a,b} satisfying 
certain restrictions. 

The following lemma gives good bounds on E{x) and F{x). 


Lemma 5.1 (Lemma 13 in [IS])- For all x > 1 we have 


( 8 ) 



n<x 


< 371-2 jQg2 ^ ^ ^ ^ ^ ^ 3.253X-^/^ 


< 67r-2xlogx + 0.787x +8.14x2/3- 0.3762. 

n<x 


One can show, using Perron’s formula and calculating residues, that 

6 6 

E{x) -^xlogx + — ( 7 r 2(27 - 1) - 12C'(2)) x, 

F(x) ~ -^xlogx + ^ (7127 - 6C'(2)) X, 

where 

— (7r2(27 - 1) - 12C'(2)) = 0.78687..., — ( 712(27 - 1) - 12C'(2)) = 1.39479 .... 

This shows that up to three decimal places, the bounds in Lemma 15.11 agree with the 
asymptotic expansions to the hrst two terms. 

Similarly, one can calculate the asymptotic order of G(x). For 3fJ(s) > 1 we have 



We can remove the contribution of p = 2 to yield 


E 

n=l,n odd 


2^uj{n) 




1 - 2 -" C^(s) 

1 + 2-C(2s)’ 
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whence, by Perron’s formula we have 


(9) 


Gix) = Y. 


n<2x,n odd 


n 


\i + 2-(*+PC(2(s+1)) s 

^{ 7 r 2(67 + 71og2) - 36C'(2)} logo:. 


vr^ 


Stt^ 


Were one to use this in the proof of Lemma 4.4 in |1] one would have 


N 


( 10 ) 


IN 


y] dnin^ - 1) < — log ^H + ... 

TT^ 


n=2 


It would take considerable effort to furnish an explicit version of fllUl) . Even if the lower 
order terms in fIlOp were negligible, one would only save a factor of two for the values of 
H and N used in our calculations. We have therefore not pursued thisS 

We content ourselves with the bound for G{x) as given in Lemma 4.3 in [4], namely 

3 0 525 

Gix) < —- log^ X + 3.1227147logx + 3.56851 H—^-. 

27r^ X 

Using this and ([H]) in ([7]) we follow the proof of Lemma 4.4 in |3] to prove 


Lemma 5.2. Let dnin) denote the number of positive integers e such that e\n and e < H. 
Then, for any N > 2 and H >1 we have 

y duin^ - 1) < iV (4 log" H + 11.1468 log77 - 0.957 + 

n=2 '' 



We now proceed to examine the number of quintuples that could arise from each of the 
triples (A)-(D). 


5.1. Case (A). This is the most damaging case in our considerations. We have r < 
((7/16)^/'^, whence, by Theorem |2] we have r < 4.611 ■ 10^® = Ra- Using Lemma 15.21 we 
hnd that the number of doubles is at most 

. JiA 

^^dRjr^-l)< 4.080-10^^. 

r=3 

Since b < {d/20y/^ < 1.9011 • 10^3 we hnd that b could have as many as 23 distinct prime 
factors. We hnd that the number of quintuples is therefore bounded by 

(11) 3 ■ 4 • 2^^ ■ 4.08 • 10^® < 8.215 • 10^7 

Since the number of possible quintuples originating from case (A) is by far the largest, 
we devote ^ to reducing this number slightly. 


^We note that Dudek has recently shown that — 1) ^ Gtt ^A^log^ N. 
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5.2. Case (B). Since 6 > 4a we have b > 2r, whence c > 4r+a. Since d > Aabc this shows 
that d > 4(r^ —l)(4r+2) > 16r^. From Theorem[2]we therefore have r < 4.147-10^® = Rb- 
By Lemma [5.21 the number of doubles {a, 6 } is at most 

.. Rb 

r=3 

Since there are at most four ways of extending a quadruple to a quintuple we find that 
the total number of quintuples is bounded above by 

( 12 ) 2.0 ■ 10^1 


5.3. Case (C). We proceed as in case 2(iii) in [16]. We consider the cases a > rj and 
a < T] and optimise over rj. In the former case, we have d > Aabc > Ar]b^^‘^ so that 
b < {d/{Ar])y^^ := N^a- Hence, by Lemma 3.3 in [5], the number of quintuples is at most 

(13) ^(logiV 3 a + 2)='-8-5-4. 

6 

When a < rj, we have b < ((i/(4a))^/® so that r‘^ = ab + 1 < a{d/{Aa)y^^ + 1. Thus 


r < 



Nsh. 


We apply Lemma [52] with H = rj and N = N^f,. Since b < (d/d)^/^ < 2.35 • 10^^ we have 
u){b) <17. Following the proof in [9] we deduce that the number of quintuples is at most 

/ 9 24 loen 44 14\ 

(14) 4-2^^ - 5 - d-Wb ^log^ 17 + 11.1468log7]- 0.957+-^ . 

\7l^ 71^7] 7] J 

We find that we can minimise the maximum of f[T3|l and f[TT|l at = 6.76 ■ 10^°. Hence 
the number of quintuples is at most 

(15) 2.41 ■ 10^^ 


5.4. Case (D). We have b < (4d/9)^/^ so that, by Theorem[21 we have b < 1.05 ■ 10^^ = 
Rd- The number of doubles {a, b} is therefore bounded by 2 We use this and 

Lemma 15.11 to prove that the number of quintuples is at most 

(16) 2.07-10^^ 


6 . Improvements to case (A) 

Here we investigate two methods. The first, in ^6.11 reduces the bound on u}{b) from 23 
to 22 , thereby saving a factor of 2 in the estimate recorded in ([TT]) . The second, in 1 16.21 
splits up the sum over b with u{b) held constant. This saves a factor of about 3.58. 
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6.1. Removing one prime factor from b. Let (pn)neN denote the sequence of prime 
numbers, and consider those b satisfying 

23 

(17) 6 o := ~ ^ uj{b)= 23. 

i=l 

We aim at enumerating all such b in flTTI) . We shall show that none of these values of b can 
appear as the second-smallest element of a quintuple. This then shows that u{b) < 22, 
and leads immediately to a saving of a factor of 2 in flTTl) . 

Suppose {a, b, c, d, e} is a quintuple. In case (A), Theorem [2] gives the bound d < UD := 
1Q®7-859_ wiien b is restricted as in flTTD we hnd that 2 divides b, since, if not, the smallest 
b can be is Y^^iPi/2 ■ p 24 > 1-18 ■ 10^"^. Continuing in this way we hnd that 2, 3, 5, 7,11 
must all divide b. 

From 4a(4a -|- 1)6^ < UD it then follows a < 7. Moreover, as the corresponding r is 
odd, ab is a multiple of 8, whence 6 = 0 (mod 8) for odd a and 6 = 0 (mod 4) for a = 2 
(mod 4). Hence, each such 6 is obtained from 6i = 6i(a) by replacing v of its factors pq, 
..., P 23 by other v primes ..., where 24 < fci < • • • < and then multiplying 
by some positive integer q such that the result is at most 

UB = UB{a, UD) := UD^^^{lQa'^ + 4 a)- 7 / 2 . 

Here 61 (a) = 46o if a is odd, 61 (a) = 26o if a = 2, 6 , and 61 (a) = 69 otherwise. 

We now present a detailed exposition of the idea sketched above. All computations 
have been performed with GP scripts [T5j. Clearly, the maximal v is determined from the 
condition 

P 24 P 25 ■ ■ ■P 23 +V ^ UB 
P 23 P 22 • • ■ P 24 -V bi 

A short computer search gives v = 3 for a = 2 or 4; n = 2 for a = 1; n = 0 for the other 
values a <7. 

Next for each m = 1, 2, ..., n we look for the largest index K = K{u) satisfying 

P 24 P 25 ■ ■ • P 22 +UPK ^ UB 
P 23 P 22 ■ ■ ■ P 24 -U 61 

and the smallest J verifying 

P24P2b • • ■P 23 +U ^ UB 

P23P22- ■ ■P25-uPj 61 

After that we determined all integers 24 < /ci < • • • < < 77 and 23 > ji > • • • > > J 

such that 

PkiPk2 ■■■Pk^ ^ UB 
PjiPj2---Pju U ■ 

Each such tuple {ki,..., ku, ji, ■ ■ ■ ,ju) gives rise to 

UBPnPj2---Pju 
_ biPkiPki ■ ■ ■ Pku _ 

candidates for the largest entry in a Diophantine couple {a, 6 }. 
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Since the bound UD = found in case (A) entails UB{a, UD) < + 

4a)for 1 < a < 7 one has 

UB{a,UD) UB{4,UD) 
bi{a) ~ 6i(4) 

Therefore, the multiplier q mentioned above must be equal to 1. 

For each value of b identihed using the above method, we are able to show easily that 
there is no corresponding quadruple. This shows that u{b) < 22. In theory there is 
nothing stopping us from playing this trick again. However, when we search for u{b) = 22 
we hnd that we could have over four thousand primes dividing b. This appears to be 
orders of magnitude harder than the u{b) = 23 case. 


6.2. Bounding b in different ranges. We have a 6 + 1 = r^. Note hrst that (i(r^ — 1) 
is even (it is odd if and only if — 1 = which implies that (r + s)(r — s) = 1 — 
a contradiction). Since d(r^ — 1 ) counts the number of divisors of — 1 it follows that 
—1) counts the number of pairs of divisors {a, b} with a < b. Now each a corresponds 
with exactly one b (and hence one pair corresponds with exactly one value of a): therefore 
|(i(r^ — 1) is actually counting the divisors a. Furthermore, note that 

(18) — 1 = a 6 > a^. 

Therefore |(i(r^ — 1) is actually counting all those a with a < \/r^ — 1. Whence for a 
hxed r we wish to count 

- 1 ). 

If r < i? then summing over r shows that the number of pairs {a, b} is at most 


(19) 


R R R 


Now, we can make a slight improvement on flT^ . Since, for case (A) quadruples we 
have b > 4a, we can improve on fllSp to show that — 1 = a 6 > 4a^. Therefore, we 
amend fll9p to show that the total number of pairs is at most 


^5^c?R/2(r'-l)- 

r=3 

One can go further than this. Let N{a,f3) be the number of quintuples with aa < b < 
(3a, for some (3 > a > 4. It then follows that for integers rrii satisfying 4 = tuq < mi < 
... < mk the total number of quintuples is bounded above by 


77(4, mi) + N{mi, m 2 ) H-h 77(mfc_i, m*,) + N{mk, 00 ), 

where N{mk,oo) means all those pairs {a, 6 } such that b > m^a. With the exception of 
N{mk, 00 ), each number is of the form 77(mj, m^+i). 
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Take rrija < b < rrij+ia. Since d > 4a6(4a6 + a + b) > 16a^6^ > 166^/(mj+i)^ we have 


( 20 ) 

We also have 


b < 


(21) — 1 = a 6 > a < R/y/mj- 

By taking rrij large we ensure that the bound on a in fl?I|) is small. We now look at u:{b) 
for b satisfying fl2U]) . We want to choose m^+i to be as large as possible such that we do 
not increase u{b). For example, when j = 0 we are considering 4a < 6 < mia. We hnd, 
using d < 7.228 x 10 ®^, that cj(b) < 14 provided that mi < 177. Also, for m 2 we hnd 
that we can take m 2 < 499686 and still ensure that cj(b) < 15. We continue in this way, 
contenting ourselves with estimates on mj that are accurate to one decimal place. We 
hnd, using Mathematica na, that we may take 

(m 3 , m 4 , ms, mg, mj, mg) = (1.7 ■ 10^ 6.4 • 10^^ 2.9 • 10^^ 1.4 • 10^°, 7.8 ■ lO^^, 4.8 • 10^^). 

We know, from 1 16.11 that there are at most 22 distinct prime factors of b. Therefore we 
have that the number of quintuples is at most 

( Ft R R \ 

2 ^^ ^ “ 1) + 2^® ^ — 1) + ■ ■ ■ + 2^^ ^ '^_r/V4.8-io27(’"^ ~ ' 

r=3 r=3 r=3 / 

We hnd that the above is no more than 


( 22 ) 1 . 177 - 102 ^. 

Using ([12]), ffTHjl . flT6|l and fl2^ we complete the proof of Theorem [TJ 
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